Low-level resource partitioning supports coexistence among functionally redundant bacteria during successional dynamics.

Members of microbial communities can substantially overlap in substrate use. However, what enables functionally redundant microorganisms to coassemble or even stably coexist remains poorly understood. Here, we show that during unstable successional dynamics on complex, natural organic matter, functionally redundant bacteria can coexist by partitioning low-concentration substrates even though they compete for one simple, dominant substrate. We allowed ocean microbial communities to self-assemble on leachates of the brown seaweed Fucus vesiculosus and then analyzed the competition among 10 taxonomically diverse isolates representing two distinct stages of the succession. All, but two isolates, exhibited an average of 90% ± 6% pairwise overlap in resource use, and functional redundancy of isolates from the same assembly stage was higher than that from between assembly stages, leading us to construct a simpler four-isolate community with two isolates from each of the early and late stages. We found that, although the short-term dynamics of the four-isolate communities in F. vesiculosus leachate was dependent on initial isolate ratios, in the long term, the four isolates stably coexist in F. vesiculosus leachate, albeit with some strains at low abundance. We therefore explored the potential for nonredundant substrate use by genomic content analysis and RNA expression patterns. This analysis revealed that the four isolates mainly differed in peripheral metabolic pathways, such as the ability to degrade pyrimidine, leucine, and tyrosine, as well as aromatic substrates. These results highlight the importance of fine-scale differences in metabolic strategies for supporting the frequently observed coexistence of large numbers of rare organisms in natural microbiomes.

Metabolite diversity among representatives of divergent Prochlorococcus ecotypes.

IMPORTANCE: Approximately half of the annual carbon fixation on Earth occurs in the surface ocean through the photosynthetic activities of phytoplankton such as the ubiquitous picocyanobacterium Prochlorococcus. Ecologically distinct subpopulations (or ecotypes) of Prochlorococcus are central conduits of organic substrates into the ocean microbiome, thus playing important roles in surface ocean production. We measured the chemical profile of three cultured ecotype strains, observing striking differences among them that have implications for the likely chemical impact of Prochlorococcus subpopulations on their surroundings in the wild. Subpopulations differ in abundance along gradients of temperature, light, and nutrient concentrations, suggesting that these chemical differences could affect carbon cycling in different ocean strata and should be considered in models of Prochlorococcus physiology and marine carbon dynamics.

Microorganisms and dissolved metabolites distinguish Florida's Coral Reef habitats.

As coral reef ecosystems experience unprecedented change, effective monitoring of reef features supports management, conservation, and intervention efforts. Omic techniques show promise in quantifying key components of reef ecosystems including dissolved metabolites and microorganisms that may serve as invisible sensors for reef ecosystem dynamics. Dissolved metabolites are released by reef organisms and transferred among microorganisms, acting as chemical currencies and contributing to nutrient cycling and signaling on reefs. Here, we applied four omic techniques (taxonomic microbiome via amplicon sequencing, functional microbiome via shotgun metagenomics, targeted metabolomics, and untargeted metabolomics) to waters overlying Florida's Coral Reef, as well as microbiome profiling on individual coral colonies from these reefs to understand how microbes and dissolved metabolites reflect biogeographical, benthic, and nutrient properties of this 500-km barrier reef. We show that the microbial and metabolite omic approaches each differentiated reef habitats based on geographic zone. Further, seawater microbiome profiling and targeted metabolomics were significantly related to more reef habitat characteristics, such as amount of hard and soft coral, compared to metagenomic sequencing and untargeted metabolomics. Across five coral species, microbiomes were also significantly related to reef zone, followed by species and disease status, suggesting that the geographic water circulation patterns in Florida also impact the microbiomes of reef builders. A combination of differential abundance and indicator species analyses revealed metabolite and microbial signatures of specific reef zones, which demonstrates the utility of these techniques to provide new insights into reef microbial and metabolite features that reflect broader ecosystem processes.

Chitin utilization by marine picocyanobacteria and the evolution of a planktonic lifestyle.

Marine picocyanobacteria Prochlorococcus and Synechococcus, the most abundant photosynthetic cells in the oceans, are generally thought to have a primarily single-celled and free-living lifestyle. However, while studying the ability of picocyanobacteria to supplement photosynthetic carbon fixation with the use of exogenous organic carbon, we found the widespread occurrence of genes for breaking down chitin, an abundant source of organic carbon that exists primarily as particles. We show that cells that encode a chitin degradation pathway display chitin degradation activity, attach to chitin particles, and show enhanced growth under low light conditions when exposed to chitosan, a partially deacetylated soluble form of chitin. Marine chitin is largely derived from arthropods, which underwent major diversifications 520 to 535 Mya, close to when marine picocyanobacteria are inferred to have appeared in the ocean. Phylogenetic analyses confirm that the chitin utilization trait was acquired at the root of marine picocyanobacteria. Together this leads us to postulate that attachment to chitin particles allowed benthic cyanobacteria to emulate their mat-based lifestyle in the water column, initiating their expansion into the open ocean, seeding the rise of modern marine ecosystems. Subsequently, transitioning to a constitutive planktonic life without chitin associations led to cellular and genomic streamlining along a major early branch within Prochlorococcus. Our work highlights how the emergence of associations between organisms from different trophic levels, and their coevolution, creates opportunities for colonizing new environments. In this view, the rise of ecological complexity and the expansion of the biosphere are deeply intertwined processes.

Advances in Chemical Analysis of Oil Spills Since the Deepwater Horizon Disaster.

Analytical techniques for chemical analysis of oil, oil photochemical and biological transformation products, and dispersants and their biodegradation products benefited significantly from research following the 2010 Deepwater Horizon (DWH) disaster. Crude oil and weathered-oil matrix reference materials were developed based on the Macondo well oil and characterized for polycyclic aromatic hydrocarbons, hopanes, and steranes for use to assure and improve the quality of analytical measurements in oil spill research. Advanced gas chromatography (GC) techniques such as comprehensive two-dimensional GC (GC × GC), pyrolysis GC with mass spectrometry (MS), and GC with tandem MS (GC-MS/MS) provide a greater understanding at the molecular level of composition and complexity of oil and weathering changes. The capabilities of high-resolution MS (HRMS) were utilized to extend the analytical characterization window beyond conventional GC-based methods to include polar and high molecular mass components (>400 Da) and to provide new opportunities for discovery, characterization, and investigation of photooxidation and biotransformation products. Novel separation approaches to reduce the complexity of the oil and weathered oil prior to high-resolution MS and advanced fluorescence spectrometry have increased the information available on spilled oil and transformation products. HRMS methods were developed to achieve the required precision and sensitivity for detection of dispersants and to provide molecular-level characterization of the complex surfactants. Overall, research funding following the DWH oil spill significantly advanced and expanded the use of analytical techniques for chemical analysis to support petroleum and dispersant characterization and investigations of fate and effects of not only the DWH oil spill but future spills.

Microbial metabolites in the marine carbon cycle.

One-quarter of photosynthesis-derived carbon on Earth rapidly cycles through a set of short-lived seawater metabolites that are generated from the activities of marine phytoplankton, bacteria, grazers and viruses. Here we discuss the sources of microbial metabolites in the surface ocean, their roles in ecology and biogeochemistry, and approaches that can be used to analyse them from chemistry, biology, modelling and data science. Although microbial-derived metabolites account for only a minor fraction of the total reservoir of marine dissolved organic carbon, their flux and fate underpins the central role of the ocean in sustaining life on Earth.

Linkages Among Dissolved Organic Matter Export, Dissolved Metabolites, and Associated Microbial Community Structure Response in the Northwestern Sargasso Sea on a Seasonal Scale.

Deep convective mixing of dissolved and suspended organic matter from the surface to depth can represent an important export pathway of the biological carbon pump. The seasonally oligotrophic Sargasso Sea experiences annual winter convective mixing to as deep as 300 m, providing a unique model system to examine dissolved organic matter (DOM) export and its subsequent compositional transformation by microbial oxidation. We analyzed biogeochemical and microbial parameters collected from the northwestern Sargasso Sea, including bulk dissolved organic carbon (DOC), total dissolved amino acids (TDAA), dissolved metabolites, bacterial abundance and production, and bacterial community structure, to assess the fate and compositional transformation of DOM by microbes on a seasonal time-scale in 2016-2017. DOM dynamics at the Bermuda Atlantic Time-series Study site followed a general annual trend of DOC accumulation in the surface during stratified periods followed by downward flux during winter convective mixing. Changes in the amino acid concentrations and compositions provide useful indices of diagenetic alteration of DOM. TDAA concentrations and degradation indices increased in the mesopelagic zone during mixing, indicating the export of a relatively less diagenetically altered (i.e., more labile) DOM. During periods of deep mixing, a unique subset of dissolved metabolites, such as amino acids, vitamins, and benzoic acids, was produced or lost. DOM export and compositional change were accompanied by mesopelagic bacterial growth and response of specific bacterial lineages in the SAR11, SAR202, and SAR86 clades, Acidimicrobiales, and Flavobacteria, during and shortly following deep mixing. Complementary DOM biogeochemistry and microbial measurements revealed seasonal changes in DOM composition and diagenetic state, highlighting microbial alteration of the quantity and quality of DOM in the ocean.

Benthic exometabolites and their ecological significance on threatened Caribbean coral reefs.

Benthic organisms are the architectural framework supporting coral reef ecosystems, but their community composition has recently shifted on many reefs. Little is known about the metabolites released from these benthic organisms and how compositional shifts may influence other reef life, including prolific microorganisms. To investigate the metabolite composition of benthic exudates and their ecological significance for reef microbial communities, we harvested exudates from six species of Caribbean benthic organisms including stony corals, octocorals, and an invasive encrusting alga, and subjected these exudates to untargeted and targeted metabolomics approaches using liquid chromatography-mass spectrometry. Incubations with reef seawater microorganisms were conducted to monitor changes in microbial abundances and community composition using 16 S rRNA gene sequencing in relation to exudate source and three specific metabolites. Exudates were enriched in amino acids, nucleosides, vitamins, and indole-based metabolites, showing that benthic organisms contribute labile organic matter to reefs. Furthermore, exudate compositions were species-specific, and riboflavin and pantothenic acid emerged as significant coral-produced metabolites, while caffeine emerged as a significant invasive algal-produced metabolite. Microbial abundances and individual microbial taxa responded differently to exudates from stony corals and octocorals, demonstrating that exudate mixtures released from different coral species select for specific bacteria. In contrast, microbial communities did not respond to individual additions of riboflavin, pantothenic acid, or caffeine. This work indicates that recent shifts in benthic organisms alter exudate composition and likely impact microbial communities on coral reefs.

Analysis of 39 drugs and metabolites, including 8 glucuronide conjugates, in an upstream wastewater network via HPLC-MS/MS.

Pharmaceutical compounds ingested by humans are metabolized and excreted in urine and feces. These metabolites can be quantified in wastewater networks using wastewater-based epidemiology (WBE) methods. Standard WBE methods focus on samples collected at wastewater treatment plants (WWTPs). However, these methods do not capture more labile classes of metabolites such as glucuronide conjugates, products of the major phase II metabolic pathway for drug elimination. By shifting sample collection more upstream, these unambiguous markers of human exposure are captured before hydrolysis in the wastewater network. In this paper, we present an HPLC-MS/MS method that quantifies 8 glucuronide conjugates in addition to 31 parent and other metabolites of prescription and synthetic opioids, overdose treatment drugs, illicit drugs, and population markers. Calibration curves for all analytes are linear (r2 > 0.98), except THC (r2 = 0.97), and in the targeted range (0.1-1,000 ng mL-1) with lower limits of quantification (S/N = 9) ranging from 0.098 to 48.75 ng mL-1. This method is fast with an injection-to-injection time of 7.5 min. We demonstrate the application of the method to five wastewater samples collected from a manhole in a city in eastern Massachusetts. Collected wastewater samples were filtered and extracted via solid-phase extraction (SPE). The SPE cartridges are eluted and concentrated in the laboratory via nitrogen-drying. The method and case study presented here demonstrate the potential and application of expanding WBE to monitoring labile metabolites in upstream wastewater networks.

Quantification of Amine- and Alcohol-Containing Metabolites in Saline Samples Using Pre-extraction Benzoyl Chloride Derivatization and Ultrahigh Performance Liquid Chromatography Tandem Mass Spectrometry (UHPLC MS/MS).

Dissolved metabolites serve as nutrition, energy, and chemical signals for microbial systems. However, the full scope and magnitude of these processes in marine systems are unknown, largely due to insufficient methods, including poor extraction of small, polar compounds using common solid-phase extraction resins. Here, we utilized pre-extraction derivatization and ultrahigh performance liquid chromatography electrospray ionization tandem mass spectrometry (UHPLC-ESI-MS/MS) to detect and quantify targeted dissolved metabolites in seawater and saline culture media. Metabolites were derivatized with benzoyl chloride by their primary and secondary amine and alcohol functionalities and quantified using stable isotope-labeled internal standards (SIL-ISs) produced from 13C6-labeled benzoyl chloride. We optimized derivatization, extraction, and sample preparation for field and culture samples and evaluated matrix-derived biases. We have optimized this quantitative method for 73 common metabolites, of which 50 cannot be quantified without derivatization due to low extraction efficiencies. Of the 73 metabolites, 66 were identified in either culture media or seawater and 45 of those were quantified. This derivatization method is sensitive (detection limits = pM to nM), rapid (∼5 min per sample), and high throughput.

Resource partitioning of phytoplankton metabolites that support bacterial heterotrophy.

The communities of bacteria that assemble around marine microphytoplankton are predictably dominated by Rhodobacterales, Flavobacteriales, and families within the Gammaproteobacteria. Yet whether this consistent ecological pattern reflects the result of resource-based niche partitioning or resource competition requires better knowledge of the metabolites linking microbial autotrophs and heterotrophs in the surface ocean. We characterized molecules targeted for uptake by three heterotrophic bacteria individually co-cultured with a marine diatom using two strategies that vetted the exometabolite pool for biological relevance by means of bacterial activity assays: expression of diagnostic genes and net drawdown of exometabolites, the latter detected with mass spectrometry and nuclear magnetic resonance using novel sample preparation approaches. Of the more than 36 organic molecules with evidence of bacterial uptake, 53% contained nitrogen (including nucleosides and amino acids), 11% were organic sulfur compounds (including dihydroxypropanesulfonate and dimethysulfoniopropionate), and 28% were components of polysaccharides (including chrysolaminarin, chitin, and alginate). Overlap in phytoplankton-derived metabolite use by bacteria in the absence of competition was low, and only guanosine, proline, and N-acetyl-D-glucosamine were predicted to be used by all three. Exometabolite uptake pattern points to a key role for ecological resource partitioning in the assembly marine bacterial communities transforming recent photosynthate.

Metabolomics as an Emerging Tool in the Search for Astrobiologically Relevant Biomarkers.

It is now routinely possible to sequence and recover microbial genomes from environmental samples. To the degree it is feasible to assign transcriptional and translational functions to these genomes, it should be possible, in principle, to largely understand the complete molecular inputs and outputs of a microbial community. However, gene-based tools alone are presently insufficient to describe the full suite of chemical reactions and small molecules that compose a living cell. Metabolomic tools have developed quickly and now enable rapid detection and identification of small molecules within biological and environmental samples. The convergence of these technologies will soon facilitate the detection of novel enzymatic activities, novel organisms, and potentially extraterrestrial life-forms on solar system bodies. This review explores the methodological problems and scientific opportunities facing researchers who hope to apply metabolomic methods in astrobiology-related fields, and how present challenges might be overcome.

AutoTuner: High Fidelity and Robust Parameter Selection for Metabolomics Data Processing.

Untargeted metabolomics experiments provide a snapshot of cellular metabolism but remain challenging to interpret due to the computational complexity involved in data processing and analysis. Prior to any interpretation, raw data must be processed to remove noise and to align mass-spectral peaks across samples. This step requires selection of dataset-specific parameters, as erroneous parameters can result in noise inflation. While several algorithms exist to automate parameter selection, each depends on gradient descent optimization functions. In contrast, our new parameter optimization algorithm, AutoTuner, obtains parameter estimates from raw data in a single step as opposed to many iterations. Here, we tested the accuracy and the run-time of AutoTuner in comparison to isotopologue parameter optimization (IPO), the most commonly used parameter selection tool, and compared the resulting parameters' influence on the properties of feature tables after processing. We performed a Monte Carlo experiment to test the robustness of AutoTuner parameter selection and found that AutoTuner generated similar parameter estimates from random subsets of samples. We conclude that AutoTuner is a desirable alternative to existing tools, because it is scalable, highly robust, and very fast (∼100-1000× speed improvement from other algorithms going from days to minutes). AutoTuner is freely available as an R package through BioConductor.

Intracellular Metabolites in Marine Microorganisms during an Experiment Evaluating Microbial Mortality.

Metabolomics is a tool with immense potential for providing insight into the impact of biological processes on the environment. Here, we used metabolomics methods to characterize intracellular metabolites within marine microorganisms during a manipulation experiment that was designed to test the impact of two sources of microbial mortality, protozoan grazing and viral lysis. Intracellular metabolites were analyzed with targeted and untargeted mass spectrometry methods. The treatment with reduced viral mortality showed the largest changes in metabolite concentrations, although there were organic compounds that shifted when the impact of protozoan grazers was reduced. Intracellular concentrations of guanine, phenylalanine, glutamic acid, and ectoine presented significant responses to changes in the source of mortality. Unexpectedly, variability in metabolite concentrations were not accompanied by increases in microbial abundance which indicates that marine microorganisms altered their internal organic carbon stores without changes in biomass or microbial growth. We used Weighted Correlation Network Analysis (WGCNA) to identify correlations between the targeted and untargeted mass spectrometry data. This analysis revealed multiple unknown organic compounds were correlated with compatible solutes, also called osmolytes or chemical chaperones, which emphasizes the dominant role of compatible solutes in marine microorganisms.

Pangenomics Analysis Reveals Diversification of Enzyme Families and Niche Specialization in Globally Abundant SAR202 Bacteria.

It has been hypothesized that the abundant heterotrophic ocean bacterioplankton in the SAR202 clade of the phylum Chloroflexi evolved specialized metabolisms for the oxidation of organic compounds that are resistant to microbial degradation via common metabolic pathways. Expansions of paralogous enzymes were reported and implicated in hypothetical metabolism involving monooxygenase and dioxygenase enzymes. In the proposed metabolic schemes, the paralogs serve the purpose of diversifying the range of organic molecules that cells can utilize. To further explore SAR202 evolution and metabolism, we reconstructed single amplified genomes and metagenome-assembled genomes from locations around the world that included the deepest ocean trenches. In an analysis of 122 SAR202 genomes that included seven subclades spanning SAR202 diversity, we observed additional evidence of paralog expansions that correlated with evolutionary history, as well as further evidence of metabolic specialization. Consistent with previous reports, families of flavin-dependent monooxygenases were observed mainly in the group III SAR202 genomes, and expansions of dioxygenase enzymes were prevalent in those of group VII. We found that group I SAR202 genomes encode expansions of racemases in the enolase superfamily, which we propose evolved for the degradation of compounds that resist biological oxidation because of chiral complexity. Supporting the conclusion that the paralog expansions indicate metabolic specialization, fragment recruitment and fluorescent in situ hybridization (FISH) with phylogenetic probes showed that SAR202 subclades are indigenous to different ocean depths and geographical regions. Surprisingly, some of the subclades were abundant in surface waters and contained rhodopsin genes, altering our understanding of the ecological role of SAR202 species in stratified water columns.IMPORTANCE The oceans contain an estimated 662 Pg C in the form of dissolved organic matter (DOM). Information about microbial interactions with this vast resource is limited, despite broad recognition that DOM turnover has a major impact on the global carbon cycle. To explain patterns in the genomes of marine bacteria, we propose hypothetical metabolic pathways for the oxidation of organic molecules that are resistant to oxidation via common pathways. The hypothetical schemes we propose suggest new metabolic pathways and classes of compounds that could be important for understanding the distribution of organic carbon throughout the biosphere. These genome-based schemes will remain hypothetical until evidence from experimental cell biology can be gathered to test them. Our findings also fundamentally change our understanding of the ecology of SAR202 bacteria, showing that metabolically diverse variants of these cells occupy niches spanning all depths and are not relegated to the dark ocean.

Hepatic metabolite profiling of polychlorinated biphenyl (PCB)-resistant and sensitive populations of Atlantic killifish (Fundulus heteroclitus).

Atlantic killifish inhabiting polluted sites along the east coast of the U.S. have evolved resistance to toxic effects of contaminants. One such contaminated site is the Acushnet River estuary, near New Bedford Harbor (NBH), Massachusetts, which is characterized by very high PCB concentrations in the sediments and in the tissues of resident killifish. Though killifish at this site appear to be thriving, the metabolic costs of survival in a highly contaminated environment are not well understood. In this study we compared the hepatic metabolite profiles of resistant (NBH) and sensitive populations (Scorton Creek (SC), Sandwich, MA) using a targeted metabolomics approach in which polar metabolites were extracted from adult fish livers and quantified. Our results revealed differences in the levels of several metabolites between fish from the two sites. The majority of these metabolites are associated with one-carbon metabolism, an important pathway that supports multiple physiological processes including DNA and protein methylation, nucleic acid biosynthesis and amino acid metabolism. We measured the gene expression of DNA methylation (DNA methyltransferase 1, dnmt1) and demethylation genes (Ten-Eleven Translocation (TET) genes) in the two populations, and observed lower levels of dnmt1 and higher levels of TET gene expression in the NBH livers, suggesting possible differences in DNA methylation profiles. Consistent with this, the two populations differed significantly in the levels of 5-methylcytosine and 5-hydroxymethylcytosine nucleotides. Overall, our results suggest that the unique hepatic metabolite signatures observed in NBH and SC reflect the adaptive mechanisms for survival in their respective habitats.

Targeted metabolomics reveals proline as a major osmolyte in the chemolithoautotroph Sulfurimonas denitrificans.

Chemoautotrophic bacteria belonging to the genus Sulfurimonas in the class Campylobacteria are widespread in many marine environments characterized by redox interfaces, yet little is known about their physiological adaptations to different environmental conditions. Here, we used liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) in a targeted metabolomics approach to study the adaptations of Sulfurimonas denitrificans to varying salt concentrations that are found in its natural habitat of tidal mudflats. Proline was identified as one of the most abundant internal metabolites and its concentration showed a strong positive correlation with ionic strength, suggesting that it acts as an important osmolyte in S. denitrificans. 2,3-dihydroxypropane-1-sulfonate was also positively correlated with ionic strength, indicating it might play a previously unrecognized role in osmoregulation. Furthermore, the detection of metabolites from the reductive tricarboxylic acid cycle at high internal concentrations reinforces the importance of this pathway for carbon fixation in Campylobacteria and as a hub for biosynthesis. As the first report of metabolomic data for an campylobacterial chemolithoautotroph, this study provides data that will be useful to understand the adaptations of Campylobacteria to their natural habitat at redox interfaces.

Formularity: Software for Automated Formula Assignment of Natural and Other Organic Matter from Ultrahigh-Resolution Mass Spectra.

Ultrahigh resolution mass spectrometry, such as Fourier transform ion cyclotron resonance mass spectrometry (FT ICR MS), can resolve thousands of molecular ions in complex organic matrices. A Compound Identification Algorithm (CIA) was previously developed for automated elemental formula assignment for natural organic matter (NOM). In this work, we describe software Formularity with a user-friendly interface for CIA function and newly developed search function Isotopic Pattern Algorithm (IPA). While CIA assigns elemental formulas for compounds containing C, H, O, N, S, and P, IPA is capable of assigning formulas for compounds containing other elements. We used halogenated organic compounds (HOC), a chemical class that is ubiquitous in nature as well as anthropogenic systems, as an example to demonstrate the capability of Formularity with IPA. A HOC standard mix was used to evaluate the identification confidence of IPA. Tap water and HOC spike in Suwannee River NOM were used to assess HOC identification in complex environmental samples. Strategies for reconciliation of CIA and IPA assignments were discussed. Software and sample databases with documentation are freely available.

Sponge exhalent seawater contains a unique chemical profile of dissolved organic matter.

Sponges are efficient filter feeders, removing significant portions of particulate and dissolved organic matter (POM, DOM) from the water column. While the assimilation and respiration of POM and DOM by sponges and their abundant microbial symbiont communities have received much attention, there is virtually no information on the impact of sponge holobiont metabolism on the composition of DOM at a molecular-level. We applied untargeted and targeted metabolomics techniques to characterize DOM in seawater samples prior to entering the sponge (inhalant reef water), in samples exiting the sponge (exhalent seawater), and in samples collected just outside the reef area (off reef seawater). Samples were collected from two sponge species, Ircinia campana and Spheciospongia vesparium, on a near-shore hard bottom reef in the Florida Keys. Metabolic profiles generated from untargeted metabolomics analysis indicated that many more compounds were enhanced in the exhalent samples than in the inhalant samples. Targeted metabolomics analysis revealed differences in diversity and concentration of metabolites between exhalent and off reef seawater. For example, most of the nucleosides were enriched in the exhalent seawater, while the aromatic amino acids, caffeine and the nucleoside xanthosine were elevated in the off reef water samples. Although the metabolic profile of the exhalent seawater was unique, the impact of sponge metabolism on the overall reef DOM profile was spatially limited in our study. There were also no significant differences in the metabolic profiles of exhalent water between the two sponge species, potentially indicating that there is a characteristic DOM profile in the exhalent seawater of Caribbean sponges. Additional work is needed to determine whether the impact of sponge DOM is greater in habitats with higher sponge cover and diversity. This work provides the first insight into the molecular-level impact of sponge holobiont metabolism on reef DOM and establishes a foundation for future experimental studies addressing the influence of sponge-derived DOM on chemical and ecological processes in coral reef ecosystems.